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Semileptonic decay of a heavy-light pseudoscalar to a light vector meson 
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We present the results of a calculation of semileptonic form factors of D and B mesons. The calculation 
uses nonperturbatively O(a) improved quenched lattice QCD and two values of the coupling. Results for charm 
mesons show reasonable agreement with experiment. The lattice results for B — > ptvi are compared with high q 2 
experimental results. 



1. INTRODUCTION 

Semileptonic decays of heavy-light mesons are 
of interest because they can be used to determine 
elements of the CKM matrix. The nonperturba- 
tive QCD information needed to understand a de- 
cay can be calculated in lattice QCD. In the case 
of the decays D — > K* and D s — ► <f> the branch- 
ing ratios have been accurately measured and the 
relevant CKM matrix element (V cs ) is well deter- 
mined by other methods EJ. These decays are a 
good test of lattice results. The decay B — > p de- 
pends on \V u b\, which is poorly known, so lattice 
results for this decay are especially useful. 

The weak matrix elements for the decay 
Pseudoscalar— > Vector (P — » V) can be param- 
eterized in terms of 4 form factors. Treating the 
V 1 and A 11 currents separately, 



(V, /c,?7|^4 M |P,p} = i(m P + m v )A 1 (q 2 )g^ l/ r) u 
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where k is the momentum of the V, p is the mo- 
mentum of the P, q = p — k and rf is the polar- 
ization vector. Results are presented here for Aq 



instead of A, where 
mp + my 
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(3) 



2. SIMULATION DETAILS 



This calculation used the nonperturbatively 
O(a) improved action and the quenched approx- 
imation. The coefficients used to improve and 
renormalize the currents were determined non- 
perturbatively in j2| . Two values of the coupling 
were used, j3 = 6.0 and j3 — 6.2. The lattice sizes 
were 16 3 x 32 and 24 3 x 48 respectively. 

The calculation used the method described in 
||. The 3 (or 1) form factors were obtained as 
free parameters in fits to three-point functions. 
The fits were repeated with different p, k com- 
binations, which correspond to different values of 
q 2 . All cases with \k\ — or 1 in lattice 
units were considered. However in some cases 
the three-point function was very noisy or did not 
seem to satisfy the assumption of the operators 
being well separated in time. One momentum 
combination was rejected at j3 = 6.2, several were 
rejected at j3 = 6.0. 

3. CHARMED MESONS 

The simulation used 3 light quark masses 
around strange, and 4 heavy quark masses around 
charm. The simulation form factors were ex- 
trapolated and interpolated to physical values of 
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quark mass. The values of the hopping parameter 
corresponding to massless, up/down and strange 
quarks were determined in Q (isospin symmetry 
was assumed) . The quark mass was improved us- 
ing the boosted perturbation theory value of b m . 

The form factors were extrapolated in light 
quark mass using the ansatz 



F(q 2 ,mi,m 2 ) = x + x\m\ + x 2 m 2 



(4) 



where F is a generic form factor, the m-i are light 
quark masses, and the free parameters. 

The Xi were determined separately for each form 
factor, momentum combination and heavy quark. 
The ansatz (jij) is a first-order Taylor expansion in 
mi, m 2 . To a good approximation the light vec- 
tor meson mass depends linearly on quark mass, 
so (|j) implicitly includes a first-order Taylor ex- 
pansion in my. 

The form factors were interpolated in heavy 
quark mass to charm. Different sensible choices 
of interpolation ansatz give almost exactly the 
same result, because the heavy quark masses sim- 
ulated are closely spaced around charm. The re- 
sults were interpolated in mp to the experimen- 
tally determined D or D s mass. This requires 
a determination of the lattice spacing, a. q 2 at 
physical quark masses was calculated using the 
experimental values of meson mass and a. 

The statistical errors are smallest in the case 
D s — > <j) because this only involves interpolation 
of the data. This is a good place to see if the 
results depend on (3. The results for D s — > <f> are 
shown in figure [l]. The (3 — 6.0 and (3 — 6.2 re- 
sults agree within errors, although unfortunately 
the (3 = 6.0 results are rather noisy. 

With form factors at just two values of (3 a 
reliable continuum extrapolation is not possible. 
To compare lattice with experiment the (3 = 6.2 
results were treated as continuum results and the 
(3 = 6.0 where discarded. The form factors were 
interpolated in q 2 using a simple pole fit, 



F(q 2 ) = 



F(0) 



1 - q 2 /M 2 



(5) 



The axial form factors have the kinematic con- 
straint A{q 2 ) = 0, which was enforced in the fits. 
The results for A 2 are too noisy to constrain a 
pole mass so this was assumed to be the same 



pole mass as for A\ as suggested by pole dom- 
inance models ||. The high q 2 part of A 2 has 
very little influence on decay rate, so the lattice 
prediction for decay rate is not biased by this as- 
sumption. 
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Figure 1. D s — » cf> form factors. 

The results of the fits to the lattice form fac- 
tors are shown in table [l]. Note that the results 
presented in this paper are preliminary. The re- 
sults shown use tq to set the scale, which gives 
a -1 = 2.913 GeV. To estimate systematic er- 
rors results were recalculated using m p to set 
the scale (a^ 1 = 2.54 GeV). The change in 
F(Q) was smaller than the statistical errors, in 
all cases. Typically the Mp are 10% smaller, 
in physical units when m p is used to set the 
scale. The lattice predictions for integrated de- 
cay rate are compared with experiment in table 
^. The leptons were assumed to be massless and 
the Zweig rule suppressed disconnected diagrams 
which contribute to the decay D s — > cj> were not 
included. There is reasonable agreement between 
experiment and the lattice results. The form fac- 
tor ratios V(0)/Ai(Q) and A 2 {0)/A 1 (0) have been 
measured for the decays D s — > <fi and D — > K*. 
In some cases there is a significant difference be- 
tween the lattice and experimental ratio. 
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Table 1 

Parameters from pole fits to the lattice form fac- 
tors. Pole mass squared is written Mj. where F 
is the form factor modelled by the pole. Masses 
are in lattice units. Errors here and in table || are 
statistical (68% confidence level). 



L> -> K* 



D 



V(0) 


0.85l| 


0.80!^ 


0.711^ 


A o (0) 


0.63! 2 . 


0.64±| 


0.59!^ 


Mo) 


0.631 2 


0.65t 2 


0.601?, 


A 2 (0) 
M 2 


0.62!!? 


0.67+^ 


0.61+^ 


0.53±| 


0.581^ 


0.571* 4 


Ml 


0.24^ 


0.28+^ 


0.31+1 


Ml 


0.69±ig 


0.9 + _j 


i-o±g 



Table 2 

Integrated decay rates using the results in table [l] 
and the PDG values for the relevant CKM matrix 







Experiment 


Lattice 


Df- 


-f <j>l+V t 


4.0 ±0.5 


5.0 ±0.3 


D+ - 




4.5 ±0.4 


5.5 ±0.5 


D+ - 


-» p°e + v e 


0.21 ±0.08 


0.19 ±0.02 



4. EXTRAPOLATION TOB^p 

To apply these lattice results to the decay 
B — ► p requires an extrapolation in heavy quark 
mass, after the light quark extrapolations. The 
extrapolation is done at fixed u) where 



1 9 
my — q 



2mpmy 
The ansatz used is || 

F(w,mp) = mp[a s (ro P )]" 



where N = -| for A u N 
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I otherwise. Ac- 



cording to heavy quark effective theory (0) has 
the correct infinite quark mass limit for lu fixed 
and close to 1. The x\ and X2 terms in (Q) are 
corrections due to the finite mass of the heavy 
quark. Some of the simulated momentum combi- 
nations have a slightly different oj for the different 
heavy quark masses. In these cases the form fac- 
tors were interpolated in q 2 to give constant u. 



In all cases the interpolation was a small shift in 
q 2 and form factor value. 

After the extrapolation B — > p form factors 
were obtained at q 2 between about 14 GeV 2 and 
qlax- The results can be extrapolated in q 2 to 
obtain form factors for the whole kinematically 
allowed range and a prediction of integrated de- 
cay rate. However the extrapolation results in a 
strongly model dependent result. Instead we in- 
terpolate the lattice differential decay rate using 



1 dT _ G 2 F q 2 [X(q 2 p 



192vr 3 m| 



(a + b{q 2 



(8) 



I Kb 1 2 dq 2 

where a, b are free parameters and X(q 2 ) is given 
in H . Note that (g) is an approximtion and is not 
expected to be reliable much below q 2 = 14 GeV 2 . 
The fit to the lattice results is good, and gives 
a = 38+J GeV 2 , b = 0t 2 2 (statistical errors). To 
estimate systematic errors the heavy quark ex- 
trapolation was repeated using the three heavi- 
est masses and X2 in p) fixed to 0. This gives 
a = 33I4 GeV 2 , b = The scale dependence 
of a, b is insignificant. 

Using \V u b\ — 3.5 x 10~ 3 and the first set of 
values for a, b gives the partially integrated de- 
cay rate Ar(14 GeV 2 < q 2 < q 2 max ) = 8.7±^| x 
10 7 s _1 . This agrees with the experimental result 
Ar(14 GeV 2 < q 2 < 21 GeV 2 ) = 7.1 ± 2.4 x 
10 7 *" 1 |. 
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